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ABSTRACT 

The Local Group starburst galaxy ICIO is the closest example of a blue compact galaxy. 
Here, we use optical gi imaging from CFHT/MegaCam and near infra-red JHK imaging 
from UKIRT/WFCAM to conduct a comprehensive survey of the structure of ICIO. We ex¬ 
amine the spatial distribution of its resolved young, intermediate and old stellar populations 
to large radius and low effective surface brightness levels. Akin to other dwarfs with multiple 
populations of different ages, stellar populations of decreasing average age are increasingly 
concentrated in this galaxy. We find that the young, star-bursting population, and the AGB 
population, are both offset from the geometric center of the older RGB population by a few 
hundred parsecs, implying that the younger star formation occurred significantly away from 
the center of the galaxy. The RGB population traces an extended structure that is typical of 
blue compact galaxies, with an effective radius of ~ 5.75 arcmins (~ 1.25 kpc). These mea¬ 
surements show that ICIO is much more extended than has previously been realized, and this 
blue compact galaxy is one of the most extended dwarf galaxies in the Local Group. The outer¬ 
most isophotes of this galaxy are very regular in shape and essentially circular in morphology. 
Based on this analysis, we do not find any evidence to suggest that ICIO has undergone a 
recent, significant, interaction with an unknown companion. 
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1 INTRODUCTION 


The Local Group is home to a large number of galaxies, the over¬ 
whelming number of which are dwarfs. The diversity of star forma¬ 
tion histories within this population is well established (e.g.. lMat^ 
1 19981 : IWeisz et^l20l4) and the origin of this diversity is the sub¬ 
ject of considerable effort that intersects several different astronom¬ 
ical topics. These range from stellar evolution to galaxy formation, 
including the role of feedback, and the interplay between secular 
and environmental dependencies on (dwarf) galactic scales. How¬ 
ever, of the more than 70 galaxies within a megaparsec or so of the 
Milky Way, only one is classified as a “starburst” system, ICIO. 


ICIO occupies an awkward position close to the Galactic plane 
(b = -3.3°), that has made study of its optical component notori- 
ously difficult. It has both a high foreground and i nternal reddenin g 
jMassev & Armandrofilll995l : 1Sanna et al.ll2008l : iKim et al1l2009h . 


but distance estimates have conve rged to place it in the range of 
700 - SOOkpc from the Milky Way l Sakai et al.|l999|:[Demers et al.l 


|2004 ITuIIv et alj|2006l: ISanna et alJl2008l : iKim et alJl2009l) . Being 
in the direction of Cassiopeia, this distance places it as an outly¬ 
ing member of the M31 subgroup. It is one of the most luminous 
of the Local Group dwarf galaxies (My - -15) and it has long 
been known to be vigorously forming stars; its overall star forma¬ 


tion rate is the highest in th e Local Group. Its high Ha luminosity 
(e.g.. lKennicutt et^l2008t) and an extremely high surface density 
of Wolf-Rayet st ars that is approximately a factor of 20 times larger 
than the LMC jMassev & Armandrofn 1 19951 ; lMasse 2 ^_Holme^ 
l2002h lead it to be classified as a starburst galaxy. Richer et alj200ll 
argue that it also meets the criteria to be classified as a blue compact 
dwarf (BCD). Whatever its classification, IClO’s star formation rate 
in comparison to our other nearby neighbours is unusually high. 


IShostak & Skillm^ ( ll989l) showed that ICIO contains signif¬ 
icant amounts of gas in a complex structure. At least two com¬ 
ponents are evidenced; the main rotating disk-like component and 
an outer, extended envelo pe that rotates in the opposite sense. 
IShostak & Skillm^ (Il989h note that this galaxy “is obviously in 
dynamical disarray”, and consider possible origins, including a 
merger with a companion of roughly equal HI mass. However, these 
authors also point out that such a morphology is not unique; they 
argue that Sextans A and NGC4449 both show similar structures 
with extended HI envelopes surrounding central rotating disks, and 
suggest that this is a sign of previously heated material being re¬ 
accreted into the center of the potential. It is notable that, since 
this work, NGC4449 has b een shown to be undergoing a merger 
with a l ower mass galaxy jMartmez-Delgado et al.l 1201 2l ; earlier 
work bv IHunter et al.lll998l suggests this could be the second such 
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merger fo r NGC4449). The comp lex gaseous morphology of ICIO 
prompted I Wilcots & Miller) ( 1 1998h to suggest that it is a galaxy still 
forming via the accretion of gas. 


The suggestion that ICIO is undergoing an interaction is an at¬ 
tractive solution to simultaneously explain its gas kinematics and its 
current intense star formation activity. Until recently, it was unclear 
if such a process was plausible for ICIO. However, ICIO is one of 
the few galaxies in the Local Group with a measured proper motion 
thro ugh measurements of m asers using the Very Long Baseline Ar¬ 
ray jBrunthaler et al.ir2007l) . New orbital solutions s uggest that its 
last M 31 pericentric passage was nearly 2 Gyrs ago ( Nideveret^ 
20^; the galaxy is now 250 kpc away from M31. iMcConnachie 
2012t) . This does not match well with the estimat ed age of the 
starbu rst (estimated to be of order ~ lOMyrs by iMassev et al.l 
I 2 OO 7 L although the exact age is uncertain and it could plausibly be 


older). Its closest known neighbor is possibly the newl y discovered 
Cassi opeia III/Andromeda XXXII dwarf spheroidal dMartin et al.l 
l2013h . at a distance of ~ 110 kpc and a luminosity of-12.3. Re¬ 

cent interactions with known Local Group galaxies therefore seem 
highly unlikely. Interactions of ICIO with an unknown - or now- 
destroyed - companion, is still plausible, particularly since ICIO 
lies in a region of sky where there can be no claims on the com¬ 
pleteness of the sam ple of Local Group galaxies (e.g.. llrwi i]|l994 
IWhiting et al.ll200% . 


Recently, iNidever et alJ l l2013l) have analysed Green Bank 
Telescope (GBT) HI surveys in the region of ICIO. They present 
evidence of an extended feature stretching 1.3 degrees to the north¬ 
west of ICIO, that displays a large velocity gradient and that ends in 
a high column-density clump with a coherent velocity. This feature 
seems to lead ICIO in its orbit (i.e., i t is no t likely due to ram- 
pressure stripping), and INidever et alJ ( 1201 3h argue that the most 
plausible explanation is that the HI is associated with an interac¬ 
tion or merger with another, unknown, dwarf galaxy, that presum¬ 
ably has a relatively low surface brightness. 

Such an interaction as proposed by INidever et alj ( 1201 3h has 
the potential to leave imprints not only on the gas distribution, but 
also on the stellar distribution at large radius, since it is feasible 
the interaction may have produced a tidal stream, or pertu rbed the 
outer regions of the galaxy. Indeed. [Loose & ThuanI ( 1 19861) showed 
that blue compact galaxies are often embedded wi thin a large 
extend ed, low surface brightness stellar distribution. ISanna et al.l 
hold) showed that ICIO is similarly embedded within an extended 
distribution of red giant stars, that extends to at least 20 arcmins. 
However, the two dimensional structure and morphology of this 
component is not yet measured. In light of recent claims, it is worth 
revisiting the extended structure of this intriguing galaxy. 


The structure of this paper is as follows: Section 2 pro¬ 
vides a description of the CFHT/MegaCam and UKIRT/WFCAM 
datasets that provide the basis for our analysis, and provides a broad 
overview of the features present in the associated color-magnitude 
diagrams. Section 3 details our attempts to deal with foreground 
contamination and to estimate a distance to ICIO using the near-IR 
dataset. Section 4 examines the spatial distribution of the young, 
intermediate and old stellar populations, with particular emphasis 
on the low surface brightness structure of the galaxy. Section 5 dis¬ 
cusses our findings, and Section 6 provides a summary of our re¬ 
sults. 


2 DATA ACQUISITION, PROCESSING AND 
CHARACTERIZATION 


The NIR dataset was obtained using the Wide Field Camera on the 
3.8m United Kingdom Infrared Telescope during a single observ¬ 
ing run in October 2008. These data were taken as part of program 
to obtain deep NIR imaging for every Local Group galaxy visible 
from M aunakea (P.I. Irwi n; p revious publications from this survey 
include ICioni et 5]|2008l and ISibbons et al.l 20121) . A detailed de¬ 
scription of WFCAM can be found in Casali et ^ ( 1200 ll) . Briefly, 
WFCAM comprises four non-contiguous Rockwell-Hawaii-II in¬ 
frared detector arrays (HgCdTe 2048x2048 pixels) that can be used 
to observe a contiguous “tile” of 0.75deg^. For ICIO, the exposure 
time of each tile in JHK was 200s, made up of 5 different dither po¬ 
sitions with 2x2 microstepping at each position to interleave a 0.2” 
sampled image. Each individual exposure was 10s. Four of these 
pointings were used to make a contiguous tile. The seeing was gen¬ 
erally very good, with an average IQ of [0.74”, 0.65”, 0.70”] for 
JHK. 


The WFCAM pipeline at the Institute of Astronomy in Cam¬ 
bridge was used to pr ocess the data , and f ollows the same proce¬ 
dures as described in ISibbons et al.l ilOld . As in that paper, the 
magnitudes quoted for each source are based on aperture photome¬ 
try, with zero-points calibrated against 2M ASS although they ar e 
not transformed into the 2MASS system dHodskin et alj |200^ . 
Rather, they remain on the WFCAM instrumental (Vega) system, 
with t he relevant transformation equations given in lHodekin et al.l 

( I 2 OO 9 I) . 

The optical dataset is archival gi (AB mags) data from the 
MegaCam wide field camera on the Canada- France-Hawaii Tele - 
scope, that has previously been presented in ISannaLt ai] (120101) . 
MegaCam is a mosaic of 36 (usable) individual 2048 x 4612 CCDs 
arranged in a 9 x 4 grid, resulting in a contiguous field of view of 
0.96 x 0.94 degrees with a pixel scale of 0.187”/pixel. The expo¬ 
sures for ICIO are the result of the co-addition of three dithered 
sub-exposures (3 x 700s in g, 3 x 400s in i). CFHT/MegaCam 
data are pre-processed by CFHT staff using the Elixir pipeline, 
which accomplishes the usual bias, flat, and fringe corrections 
and also determines the photometric zero point of the observa¬ 
tions. These images were then processed using a ve rsion of the 
Cambridge photometry pipeline dlrwin & LewisilLooI!) adapted for 
CFHT/MegaCam observations, and which has been used exten- 
sively for the Pan-Andro meda Archaeological Survey (PAndAS; 
IMcConnachie et ^ l2009l) . The pipeline includes re-registration, 
stacking, catalog generation, and object morphological classifica¬ 
tion, and creates band-merged gi products for use in the subsequent 
analysis. The image quality of the i-band data is consistently bet¬ 
ter than the g-band data (~ 0.8 - 0.9” instead of ~ 1.0 - 1.1”). 
We found that a deeper and more consistent stellar catalog was 
produced when we used objects identified in the i-band to force 
photometry in the g-band. The object classification fo the CFHT 
data is therefore taken from the i-band only. 

The final WFCAM merged dataset contains a total of 224 679 
sources, of which 130581 are confidently classified as stellar 
(within 2 sigma of the stellar locus) in at least 2 of the three bands. 
The final MegaCam catalog contains a total of 290 720 sources, of 
which 210 156 are confidently classified as stellar (within 2 sigma 
of the stellar locus) in the i-band. 

Figure 1 shows the derived photometric errors for stars in each 
of five bands as a function of magnitude. Individual sequences cor¬ 
responding to individual detectors are particularly noticeable for 
the MegaCam i bands. In general, our photometric errors are better 
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Figure 1. Photometric uncertainties as a function of magnitude for the two optical bands from the CFHT/MegaCam data (left panels) and the three NIR bands 
from the UKIRT/WFCAM data (right panels). 
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Figure 2. Luminosity functions of all stellar sources identified in the two optical bands from the CFFlT/MegaCam data (left panels) and the three NIR bands 
from the LfKIRT/WFCAM data (right panels). For the NIR luminosity functions, stellar classification has been conducted on each band independent of all other 
bands. For the optical bands, object identification and classification was conducted exclusively using the /-band. Vertical dashed lines indicate the turn-over 
in the luminosity functions, indicating significant incompleteness at these magnitudes. 


than 10% for [g, /, J, H, K] ^ [25.0,23.2,19.8,19.2,18.5]. Figure 2 
shows the luminosity functions for stellar objects in each of the 
5 bands to illustrate the general depth of the data in each band. 
Vertical dashed lines at [g, i, J, H, K\ ^ [26.3,23.8,20.3,19.4,18.9] 
indicate the turn-over in the luminosity functions where the incom¬ 
pleteness has become sufficiently large to dominate over the in¬ 
creasing star counts at faint magnitudes. The analyses presented 
in this paper do not require precise knowledge of the shape of the 
incompleteness function as a function of magnitude. 

Finally, we correct the observed magnitudes for reddening. 
Normally, we would do this on a star-by-star basis by cross- 
ref erencing to the E(B-V ) maps derived from IRAS/DIRBE data 
by ISchleeel et al.l ( Il998h . However, the low latitude - and conse¬ 
quently extremely high expected foreground reddening - of 1C 10, 
means t hat we are working in a regime in which the reddening 
maps of ISchlegel et alj ( 1199 Sh are known to be inaccurate. Eor ref- 
erenc e, the mean reddenin g measured for stars in our 1C 10 fields 
from ISchleeel et al.l (Il998l) is E(B-V)~ 0.6 5 mags (after applying 
the correction from [Bonifacio et alJ ( |2000|) that is recommended 
in regions of moderately high reddening). However, independent 


estimates of the total reddening (i.e., internal plus foreground) in 
1C 10 by oth er authors suggest the value is nearer 0.7 - 0.8 (e.g. 
0.77 ± 0.07 jRicher eTaDl200ll: 0 .78 ± 0.06. ISanna etaDl2008t but 
also 0.98 ±0.06, Kim et al.l2009l) . Given these uncertainties, adopt¬ 
ing a spatially varying value on a star by star basis is not warranted, 
and we adopt a con stant reddening value of E(B-V)= 0.77 ± 0.07 
([Richer et alj|200lh . In what follows, adoption of different redden¬ 
ing assumptions do not have any qualitative changes on our results. 


3 COLOR-MAGNITUDE DIAGRAMS AND 

IDENTIFICATION OF THE STELLAR POPULATIONS 

Figure 3 shows the extinction-corrected color magnitude diagrams 
of stellar sources identified in the full CFHT/MegaCam (left panel; 
di < 0.5, dg < 0.3 mags) and UKIRT/WFCAM (right panel) 
fields. Stellar isochrones corresponding to a 13 Gyr population with 
[Fe/H] = - 2.0, -1.5, and -l.Od ex from the Dartmouth stellar evolu¬ 
tion library dPotter et al.l2008l) are overlaid on both CMDs, shifted 
to a distance of 750 kpc, to highlight the approximate location of 
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Figure 3. Color-magnitude diagrams (CMDs) of stellar sources in the entire CFHT/MegaCam field of view (left panel; di < 0.5, dg < 0.3 mag) and the 
UKIRT/WFCAM field of view (right panel). These are plotted as Hess diagrams with square-root scaling to better show the structure in these diagrams. St ellar 
isochrones corresponding to a 13 Gyr population with [Fe/H] = —2.0, -1.5, and -l.Odex from the Dartmouth stellar evolution library jPotter et alj|2008l) are 
overlaid on both CMDs, shifted to the approximate distance modulus of ICIO. The overwhelming fraction of stars in both CMDs are foreground stars in the 
Galaxy, that swamp the ICIO population. 
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Figure 4. The left panel shows the optical color-magnitude diagram (CMD) from the CFHT/MegaCam for stars located in a circular annulus between 2 and 3 
arcmins from the center of ICIO. This annulus has been chosen to sample an area where the ICIO population is strong in comparison to the Galactic foreground, 
while avoiding the very central regions where crowding is significant. The right panel is a reference field, at a random location in the field, offset from ICIO 
and with the same total area. It shows the presence of the Galactic foreground, but an absence of any other stellar populations. Both CMDs show only stellar 
sources with di < 0.5, dg < 0.3 mag 
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Figure 5. Near-IR color-color diagram of stellar' sources in the UKIRT/WFCAM field with uncertainties less than 0.1 mags in each filter. Red points shows the 
location of stars located more than 20 arcmins from the center of ICIO, and will predominantly consist of Galactic foreground dwarfs Black points show the 
positions of stars within 2 arcmins of the center of ICIO, and will consist predominantly of giants in ICIO. The overwhelming majority of the candidate dwarfs 
lie blue-ward of (/ — H)q = 0.7, whereas only a small fraction of the candidate giants are blue-ward of this limit. 



red giant branch stars in these plots. Both panels of Figure 3 are 
dominated not by stars in ICIO, but by the Galactic foreground 
that lies along our line of sight to this very low-latitude galaxy 
{b = -3°). This foreground is visible in the NIR CMD as two 
“fingers” extending over the full magnitude range, centered around 
(J - K)o 0.2 and (J - K)o ^ 0.6. The bluer of these features con¬ 
sists predominantly of Galactic main sequence turn-off stars (F-K 
dwarfs) at a range of distances along our sight-line in the disk and 
halo. The redder feature consists mostly of disk K d warfs, with 
some disk K giants fe.g. INikolaev & Weinberell2000h . The same 
features are also present in the optical CMD, where the turn-off 
stars are the vertical sequence at (g - i)o =: 1.0, and the disk dwarfs 
form the red cloud of points around (g-i)o ~ 2-3. In both CMDs, 
there is a notable color gradient in these features, with the fainter 
objects appearing redder than the brighter objects. This is a distance 
effect, since the Galactic sources further from us have a higher in¬ 
tegrated reddening along our sight-line than the closer sources. 

To better increase the contrast between stellar populations in¬ 
trinsic to ICIO and the Galactic foreground, the left panel of Fig¬ 
ure 4 shows the optical CMD for stars within 4 arcmins of the cen¬ 
ter of ICIO (with di < 0.5, dg < 0.3 mags). Detections within the 
inner 3 arcmins have been removed from this CMD, since detec¬ 
tions in this region are highly unreliable due to severe crowding 
issues. The right panel of Figure 4 shows the optical CMD of stars 


in a random region of the MegaCam field that samples the same 
total area as the left panel and is well away from the expected lo¬ 
cation of ICIO. While the foreground populations that so dominate 
in Figure 3 are still present in both panels, the field centered on 
ICIO has a clear excess of stars compared to the reference field. In 
particular, a large number of stars with (g - i)o < 0.5 are present, 
and correspond to young main sequence and blue supergiant stars 
in this galaxy. Further, there is a very large number of points, with 
colors typically redder than (g - i)o >0.8 and if, > 20. These are 
mostly red giant branch stars, as well as younger, more luminous 
AGB stars. It is also clear, however, that these stars do not occupy 
a clear, concentrated locus that follows the isochrones, but rather 
the points cluster in a “cloud” around the approximate location of 
the isochrones. This is due to the uncorrected effects of extinction 
discussed earlier. For this reason, we consider it wholly unreliable 
to draw any significant conclusions from the colors or luminosi¬ 
ties of these optical data (e.g., such as inferring the TRGB or the 
metallicity via the color of the RGB). 

Considering now the near-IR data, we note that it is possible to 
remove the overwhelming majority of for eground (dwarf) stars us- 
ing a near-IR color-cut, a s demonstrated in iGullieuszik et aO ( 120081) 
and lSibbons et alj ( l2012h . Figure 5 shows a (J-H)o, (H-K)o color- 
color diagram of all stellar sources in the UKIRT/WFCAM field. 
Only sources with uncertainties less than O.lmags in each filter are 
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Figure 6. NIR CMD of stellai' sources in the UKIRT/WFCAM. The red density plot shows all sources with {J — H)o < 0.7, and which is dominated by the 
Galactic foreground population. Black points are stars with (7 - 7/)o > 0.7, and are mostly genuine gi ant stars in ICIO. S tellar isochrones corresponding to a 
13 Gyr population with [Fe/H] = —2.0, -1.5, and -1.0 dex from the Dartmouth stellar evolution library tootter et alj2008h are overlaid. Stars brighter than the 
RGB are mostly asymptotic giant branch stars (AGB) stars. The cloud of stars redward of {J — K)o ~ 1.4 between 16 < ^"0 < 18 are C-rich AGB stars, and 
0-rich AGB stars are visible in the bright vertical sequence at {J - K)o ~ 1 


plotted. Red points show the location of all stars located more than 
20 arcmins from the center of ICIO, and will predominantly con¬ 
sist of Galactic foreground dwarfs. Black points show the positions 
of all stars within 2 arcmins of the center of ICIO, and will consist 
predominantly of giants in ICIO. The overwhelming majority of the 
candidate dwarfs lie hlue-ward of (J - H)o = 0.7, whereas only a 
small fraction of the candidate giants are hlue-ward of this limit. 
We can adopt (J - H)o = 0.7 as a color-cut to remove the majority 
of the foreground population (specifically, 96% of all foreground 
stars lie hlue-ward of this limit in Figure 4). The removal of the 
small number of genuine ICIO stars hy this cut makes no differ¬ 
ence to the subsequent analysis. 

Figure 6 shows the same near-IR CMD as in the right panel 
of Figure 3. The black points indicate all stellar sources with (J - 
H)q > 0.7 and are mostly members of ICIO. The red density plot 
shows the location of all stars with (J - H)q < 0.7, and are mostly 
foreground stars, as discussed previously. This CMD now clearly 
shows distinct populations of AGB stars. In particular, O-type AGB 
stars are visible in the broad, vertical sequence at around (J- K)o = 
1, and C-type AGB stars are visible at around 16 < Alo < 18 in a 
sequence that extends red-wards in color to (J - K)q ~ 2. 

In addition to the AGB population, RGB stars are also present 
in Figure 6 and occupy broadly similar colors as the (O-type) AGB 
stars, but at fainter magnitudes. Figure 7 shows the A'p-band lu¬ 
minos ity function for the ICIO stars in Figure 6. iBellazzini et al.l 
1I2OO4I) provide a calibration for the tip of the red giant branch 
as a distance indicator using near-IR wavelengths. At these wave¬ 
lengths, it is a much stronger function of metallicity than in the 
/-band, which complicates its usage for complex stellar popula¬ 
tions where stars occupy a broader metallicity range. Nevertheless, 
identification of plausible values of the tip in the K band provides 
us with a rough estimate of its distance. 


Unlike equivalent luminosity functions for dwarf galaxies 
that do not have a significant AGB population, the discontinu¬ 
ity in counts due to the tip of the RGB is not clearly defined. 
In this specific situation, the application of edg e-detection tech¬ 
niques (e.g., the Sobel edge detection algorithm. lUee et ^1 19931 : 
ISakai et al.l 19961) or bavesian or maximum likelihood techniques 
(e.g.. lMakarov et aDl2006l : IConn et al.ll201lh . that do not explicitly 
account for the (difficult to model) AGB population, do not help 
identify any particular feature as due to the tip. By joint inspection 
of the luminosity function and CMD, however, we can identify a 
plausible range of values. The largest discontinuity in the luminos¬ 
ity function is at Al 18.4, and marks the faintest plausible loca¬ 
tion we identify as the tip; at brighter magnitudes, there is a steady, 
roughly constant rise in star counts from Kq ^ 16.0 to A'o 17.95 
(that corresponds primarily to the bright AGB populations in the 
CMD) and then a change in slope to a steeper rise in counts with 
magnitudes. This change in slope could plausibly be due to the on¬ 
set of the RGB. 

In the 2MASS filter system, IBellazzini et alj (l2004b derive 
the following relationship between TRGB magnitude and overall 
metallicity for K\ 


j^TRGB ^ _0.64[M///] ■ 


6.93 


( 1 ) 


We a dopt a metallicity value of [Fe/H ] = [M/H] = -1.3 dex for 
ICIO jTikhonov & Galazutdinovall20()^ . which yields an absolute 
K-band magnitude of -6.1 mags (the correction between the WF- 
CAM and 2M ASS /f-band systems is neg l igible ). The metallic¬ 
ity estimated bv lTikhonov & Galazutdinov^ ( I 2 OO 9 I 1 is based on the 
RGB and so is appropriate for our estimate, although we also note 
that higher metallicity values exist based on planetary nebulae and 
HII r egions, that sugges t its metallicity is closer to th at of the SMC 
(e.g.. lRicher et al.ll200ll : lMagrini & Goncalvesll2009h . A change in 
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Figure 7. band luminosity function of all stellar sources shown in Figure 6 with > 0.7. The large number of AGB stars in ICIO makes determination 

of the onset of the RGB very difficult. The likely TRGB range is indicated; the bright-end limit to this range is marked by a change in slope of the luminosity 
function, and the faint-end limit is set by a large jump in star counts. 


the metallicity of ~ 0.2 dex will change the absolute magnitude 
of the TRGB by ~ 0.1 mags. For the suspected TRGB location 
of 17.95 < K ^ 18.45, this corresponds to a distance modulus 
of 24.05 < {m - M)q < 24.55 (650 ^ D < SlOkpc). This is 
consiste nt with previous estim ates (e.g. , 660 kpc, Sakai et al. 1 19991: 
741 kpc . iDemers et al] l2004 795kt)c.lTullvetalJl2006l : 798 kpc. 
Is anna et al.ll2008l: 715 kpc.lKim et al.ll2009h. 


4 STRUCTURAL PROPERTIES 
4.1 Stellar density maps and shapes 


equivalent contour map for the UKIRT/WFCAM dataset, where 
only candidate AGB stars have been selected, using the foreground 
color-color cut described earlier and requiring Kg < 17.95 (i.e., 
brighter than the expected magnitude of the TRGB). Pixel sizes and 
smoothing is the same as for the optical dataset. Contours are set at 
2-, 5-, 10-, 20-, 40-, 60-, 80 - tr above the background level. 

We quantify the two-dimensional shape of the stellar popula¬ 
tion as traced in the contour maps of ICIO via a moments analysis 
of the density distribution. Specifically, we calculate the center of 
gravity (x^, y„), position angle, 6 (measured east from north), and 
eccentricity, e, via 


We now examine the spatial distribution of the three main stellar 
populations identified in our ICIO data, namely young blue stars, 
AGB stars (O-type and C-type), and RGB stars. AGB and RGB 
stars are present in both the optical and near-IR CMDs, but the 
bright AGB stars are most uniquely identified in the near-IR CMD. 
RGB stars, on the other hand, can be identified as being the ma¬ 
jority of stars fainter than the TRGB in the near-IR CMD, but the 
depth of the WFCAM data means that this dataset is sub-optimal 
for examining the distribution of this population. In the optical data 
it is not possible to disentangle the AGB from the RGB due to the 
confusion introduced by the reddening, but weighed against this 
consideration is the fact that the optical data is considerably deeper, 
thus identifying considerable more RGB stars than in the near-IR. 

The left panel of Figure 8 shows the spatial distribution of 
stars in the CFHT/MegaCam optical dataset. The grey-scale shows 
a density plot of candidate RGB/AGB stars, selected to have 
(g - i)o > 0 (to avoid young stars) and ig > 20.4 (to reduce fore¬ 
ground contamination). The large gaps that are present between 
the first-and-second, and third-and-fourth, rows of CCDs in Mega- 
Cam have been artificially filled in this image using estimates of 
the local background. Pixel size is 0.3 arcmins, smoothed with a 
Gaussian kemal with a dispersion of 1.5 pixels. Blue points in¬ 
dicate candidate young stars, selected to have (g - i)o < 0 and 
ig < 21. The right panels show the same density map as shown in 
the grey-scale, but now as a contour map. The lowest level con¬ 
tour is set at 2 - cr above the background, then 5 - cr, increasing 
to 35 - (T in increments of 5 - cr. No contours are shown for the 
very central regions, since here there is extreme crowding and we 
focus our analysis on the outer regions only. Figure 9 shows the 


H, Xjl (x„ y,) _ 'ZiyiljXi, yd 

~ Y.il{x„yd 

1 2cr„ 

9 = — arctan-^— 

^yy ^XX 

- O-yyf + 40 -^ 

e = - (2) 

^ XX “I" ^yy 

where /(x,, y,) is the intensity of the i’’’ pixel and cr„, cTyy and <r„ 
are the intensity-weighted second moments of the distribution. The 
eccentricity is related to the ellipticity, e = I -bja, via 


6 = 1 - 



(3) 


These quantities are calculated as a function of radius using an it¬ 
erative technique. That is, for all pixels lying within a given radius 
r of a first estimate of the center, we calculate the corresponding 
value of xg,yg,e,6. We then repeat the calculation but this time 
only using pixels lying within the appropriately positioned ellipse 
described by these parameters (where r now describes the semi¬ 
major axis of the ellipse). This process is repeated until conver¬ 
gence for each value of r. Results are shown in Figure 10 for the 
optical data (left panel) and the near-IR data (right panel). For each 
stellar population, the estimates within the inner few arcminutes are 
either unreliable or impossible due to the effects of crowding; these 
regions are marked by the shaded areas in Figure 10. In these plots, 
Axo, Ayo corresponds to the offset from an adopted center for ICIO 

of (00'’20"’17.3',59°18'14”). 
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Figure 8. Left panel: blue points indicate the positions of candidate young stars in ICIO, identified in the CFHT/MegaCam data using (g — i)o < 0 and /q < 21. 
Grey-scale shows candidate RGB/AGB stars, selected as described in the text. Pixels are 0.3 X 0.3 arcmins, convolved with a Gaussian kernal of dispersion 1.5 
pixels. The density map is shown with a square-root scaling. Right panel: contour map of the RGB distribution, where the density map is the same as in the 
right panel. Contours are set at 2-, 5—, 10—, 15 — ...35 — cr above the background level. The center of ICIO is positioned at cr = OOhlOsll.3m, S = 59° 18'14”. 



^ (arcmins) 

Figure 9. Contour map of the AGB distribution from the UKIRT/WFCAM data, selected as described in the text. Pixels are 0.3 X 0.3 arcmins, convolved with 
a Gaussian kernal of dispersion 1.5 pixels. Contours are set at 2-, 5—, 10—, 20-, 40-, 60-, 80 — rr above the background level. The center of ICIO is positioned 
at O' = 00/i20xl7.3m, 6 = 59° 18' 14”. 


Table 1 lists the average values for the various parameters for 
each stellar population. This table also includes an estimate of the 
position of the center of the young population. Due to the limited 
number of young stars, the center of the population is estimated as 
the median position of all stars within 10 arcmins of the nominal 
(literature) center of ICIO. We find that the center of the young 
stars is within arcseconds of the nominal center, which is entirely 


consistent with the expectation that the recognised “center” of IC10 
coincides with the center of the stellar population that completely 
dominates the light. The scale radius has been crudely estimated as 
the median distance of these stars from the center of the population. 
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Figure 10. Variation in the shape of ICIO as a function of semin-major axis, as estimated from the RGB/AGB distribution using the CFHT/MegaCam data (left 
panels) and the AGB distribution using th UKIRT/WFCAM data (right panels). The variation in the center of the best fitting ellipse (xo,yo; top two panels), its 
ellipticity (third panel) and position angle (measured east from north; fourth panel) are plotted as a function of the semi-major axis. The grey region in the left 
hand plot indicates the region where estimates are less reliable due to crowding effects. See text for details. 



(Ax, Ay) 

arcsecs pc (D/750kpc) 

O 

6 

1-b/a 

n 

arcmins 

Te 

pc (£)/750kpc) 

RGB/AGB (optical) 

(62.1,-22.7) 

(226,-83) 

48 

0.15 

0.75 ± 0.06 

5.75 ±0.11 

1254 ±24 

AGB (near IR) 

(13.6, 0.7) 

(49, 3) 

100 

0.07 

1.20 ±0.17 

3.22 ± 0.30 

702 ± 65 

Young (optical) 

(1.6, 1.3) 

(6, 5) 

— 

— 

— 

~ 1.75 

~ 382 


Table 1. Shape and structural parameters for each of the main stellar populations identified in ICIO. Positional offsets are measured relative to the literature 
center for ICIO of (00*20'" 17.3®, 59° 18'14”), derived from the 2MASS Extended Source catalog. 


4.2 Radial profiles 

We now calculate the radial profiles of the RGB/AGB population 
shown in Figure 8, and the AGB population shown in Figure 9. 
As shown in Figure 10 and Table 1, both populations are rela¬ 
tively circular. To facilitate comparison between the two profiles, 
we therefore adopt circular annuli centered on the nominal center 
of ICIO from the literature (00*20'" 17.3', 59°18'14”). Background 
levels are estimated in the outer parts of the fields, and subtracted 
from the profiles . We a lso apply a crowding correction, following 
lirwin&Trimbleld 19841) : 

/'-/|l+2 /A' + y/A'l..J . (4) 

/ is the observed number density of all images, /' is the actual 
number density and A' is the typical area of the image, where the 
radius of the image is closely approximated by the seeing. This cor¬ 
rection ignores second-order effects such as the actual shape of the 
luminosity function. Clearly, this correction is only applicable in 
partially crowded fields, and breaks down once crowding becomes 
sufficiently severe. 

Radial profiles for ICIO are shown in Figure 11. The left panel 
corresponds to the RGB/AGB population (optical), and the right 
panel corresponds to the AGB (near-IR) population. Error bars in¬ 
dicate the combined uncertainty of the counts and the background 
estimate. Only the data for regions beyond 2 arcminutes are shown. 


since in both datasets, crowding is severe interior to this radius. 
Both profiles are remarkably regular, with no obvious breaks or 
discontinuities. 

We fit a Sersic function to each radial profile, given by 


l{r) = /o exp 




(5) 


/o is the central intensity, is the effective radius, n is the Sersic 
index, and b„ is a quantity that depends upon n and is well ap- 
proximated by h„ ^ 1.9992n - 0.3271 for 0.5 < n < 10 (e.g., 
iGraham & Drive j l2005h . Best fit values for each profile are in¬ 
cluded in each panel, and summarised in Table 1. 


5 DISCUSSION 

5.1 The distribution of stars through time 

Over the past decade, it has become increasingly clear that 
dwarf galaxies exhibit complex spatial structures i n terms of age, 
metall icity and chemical variations. For example, iHarbeck et al.l 
(l200lh showed that population gradi ents are common in the Lo¬ 
cal Group dwarf galaxy population. iTolstov et al.l (|2004|) subse¬ 
quently showed that the prototypical Sculptor dwarf spheroidal 
galaxy consists of two dynamically and chemically distinct com¬ 
ponents, wherein the younger, more metal-rich population is more 
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radius (arcmins) radius (arcmins) 

Figure 11. Radial profiles of the RGB/AGB population from the CFHT/MegaCam data (left panel) and of the AGB population from the UKIRT/WFCAM 
data (right panel), along with best-fitting Sersic profiles (dashed lines). For both panels, the central regions where crowding is significant have been omitted 
from the fits. 


centrally concentrated than the older, more metal-poor, component. 
It is now generally accepted that different stellar populations within 
the same dwarf galaxy can have d i stinct spatial distribu t ions (e.g., 
Battaglia et alj I 2 OO 6 I . I 2 OIIL1201 2 l : [Bernard et al.l I 2 OO 8 I : iLee et al] 
20091 among many others). 

Based on our analysis, we confirm that ICIO exhibits the be¬ 
haviour of stellar populations seen in many other Local Group 
gaaxies. Firstly, there is a clear age gradient, such that the young 
main sequence and blue loop stars are the most concentrated of 
the stellar populations identified in our data, and lik ely have ages 
in the range of tens to a few hundred Myrs (see also lMassev et al.l 
l2007h . Secondly, bright AGB stars, that may have ages as young 
as ^ 1 Gyr, have an effective radius that is of order twice as ex¬ 
tended as the youngest stars. Finally, the RGB stars, that are gener¬ 
ally old (>> 2-3 Gyrs), have an effective radius that is roughly 4 
times as extended as the youngest population. While typical of Lo¬ 
cal Group galaxies, this behaviour is also typical of other systems 
classified as blue compact galaxies, that are known to be embedded 
within extended low s urface brightness distributions of old stars 
jLoose & Thuanl[l986h . 

While each stellar population has a different characteristic 
scale radius, the center of the young population in ICIO is clearly 
offset from the oldest RGB population. Unsurprisingly, the centroid 
of the young stars is very close to the literature value for the central 
coordinates of ICIO, since the latter are based on the centroid of 
the light distribution (measured from the 2MASS All-Sky Release 
Point and Extended Source Catalogs and Image Atla^, which is 
dominated by the contribution of the starburst population. The cen- 


* http://www.ipac.caltech.edu/2mass/releases/allsky/index_extmis.html 


troid of the AGB population is also close to this position. However, 
the center of the oldest stellar population - the RGB - is very sig¬ 
nificantly different from these coordinates. As listed in Table 1, the 
offset of the RGB centroid to that of the young stars is ~ 240 pc, 
for an adopted ICIO distance of 750kpc. To put the significance of 
this number in context, the half-light radius of the Sculptor dwarf 
spheroidal galaxy is ~ 280 pc. 


As the oldest stellar population, RGB stars are expected to 
be the most relaxed of the stellar populations that we examine. As 
such, their centroid is likely close to the centroid of the net gravita¬ 
tional potential of the galaxy. Star formation in dwarf galaxies need 
not be a uniform and symmetric process. Indeed, an irregular dis¬ 
tribution of young stars is a defi ning morphologic al characteristic 


of Magellanic irregular galaxies I Holmberall958l). and IC10 wa s 
first classified as such a system by de Vaucouleurs & Abies! (Il965h . 


iBastian et alJ l l201 ih show, through tracing blue loop stars, that star 
formation in dwarf galaxies is highly substructured on scales of 
tens to hundreds of parsecs, and that such subs tructures can persist 
for po tentially hundreds of megayears (see also lDohm-Palmer et al.l 
1 2002 !) . In this picture, it is not surprising that the nearest galaxy 
with the highest star formation rate has ongoing star formation that 
is not centered in the gravitational potential of the galaxy. Indeed, 
if the star formation in ICIO was induced through an interaction, 
then it is certainly not clear that the star formation will be centered 
with respect to the net gravitational potential. 


Finally, we suggest that it is important that the difference be¬ 
tween the center of the young stellar population in ICIO and that 
of the old population is considered when conducting dynamical 
analyses of ICIO in which a model of the gravitational potential 
is being used. Specifically, large systematic uncertainties could be 
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introduced if models fail to properly account for the fact that the 
center of mass of this galaxy is not the same as its center of light. 

5.2 The outer structure and interaction history of ICIO 

iLoose & Thuanl ( Il986h demonstrated that blue compact galaxies 
consist of at least two components; the central, star-bursting pop¬ 
ulation and an old, extended, l ow surface brightness surrounding 
population. ISanna et^ l l201(]h showed that red giants surround 
ICIO extending to large radius, and here we measure the shape 
parameters describing this stellar component. In addition, we also 
show that stars with ages intermediate between the youngest and 
oldest populations have a spatial extent that is also intermediate 
between these populations. 

The spatial extent of the outer component of RGB stars in 
ICIO is worthy of comment. Firstly, the best-fit Sersic profile to 
the RGB distribution yields n = 0.75 ± 0.06. In this respect, 
this component of ICIO is typical of many dwarf galaxies for 
which the stell ar distribution falls o ff more slowly than an ex¬ 
ponential fe.g.. iFerrarese et n]|2006h . With an effective radius of 
Tg = 5.75 ± O.llarcmins (~ 1.2kpc for an assumed ICIO dis¬ 
tance of 750kpc), very few dwarf galaxies in the Local Group are 
as large as this old compone nt of ICIO. Comparis on with the up¬ 
dated data tables presented in lMcConnachiel ll2012H^ show that (ex¬ 
cluding the LMC, SMC and M33) there are only 6 Local Group 
galaxies with a robust measurement of their size that places them 
with r;, > Ikpc (measured on the major axis). Of these, the Sagittar- 
ius dwarf spheroidal is known to be exc eptionally tidally perturbed 
dlbata et alj200ll ; lMaiewski et al.l2003l and many others). Four oth¬ 
ers in this sample are M31 satellites (Andromeda II, XIX, XXIII 
and XXXII), and each of these is also a dwarf spheroidal. ICIO is 
also a satellite of M31, but is not classed as a dwarf spheroidal due 
to its young stars and gas. The remaining galaxy in this sample is 
WLM, a relatively isolated dwarf irregular galaxy of comparable 
luminosity to ICIO. 

It is important to note that, in addition to a large spatial ex¬ 
tent, the RGB population (and also the AGB population traced with 
the near IR data) has a very regular morphology. This is visible 
in both panels of Figure 8, where the isophotes do not show any 
obvious signs of distortion, or any features that might suggest the 
presence of stellar streams or shells. Figure 10 shows the “intensity- 
weighted” shape of the galaxy as a function of radius as traced by 
the stellar populations. At large radius, ICIO is well approximated 
as circular with no evidence of significant isophote twisting or other 
irregular behavior. There is similarly a lack of any notable features 
at large radius visible in the radial profiles in Figure 11. In both the 
optical and near-IR data, the decline in star counts matches well 
the Sersic profile fits. There is no sign that the profiles truncate 
at any point, and so it is probable that stars will be found at even 
larger radius than probed by our data. There is also no clear sign of 
any “extra-tidal stars”, that is, stars located at large radius in num¬ 
bers exceeding that expected from a simple extrapolation from the 
numbers at smaller radii. The presence of such stars has been sug¬ 
geste d in the past to be due to tidal effects dlrwin & Hatzidimitrioul 
Il995h . although other explanations such as multiple components 
are also possible. We note that there is a slight deviation of the 
RGB star counts from the fitted profile in the left panel of Fig¬ 
ure 11 at around 15 arcmins; this could plausibly be interpreted as 
marking the truncation of ICIO with the stars at larger radii being 

^ Available at http://www.astro.uvic.ca/~alan/Nearby_Dwarfs_Database.html 


“extra-tidal”. However, such an interpretation is at the limits of the 
data, and we consider it unlikely given that no similar feature is ob¬ 
served in the near IR data. We conclude that, given the data in hand, 
the outskirts of ICIO are unremarkable in their regularity, which is 
itsel f of interest in c ompar ison to the HI distribution. 

IPemers et al.l l l2004l) analyse the structure of ICIO using 
CFHT/CFHT12K wide field photometry (42' x 28') in the Mould 
I and R filters and narrow-band CN and TiO filters, that they use 
to investigate the spatial distribution of Carbon stars and the rest 
of the red giant population. They conclude that the C stars follow 
the bulk of the rest of the stellar population of ICIO, and that both 
are well described by power law profiles. Our Sersic profiles in the 
radial range o f 5-14 arcmins from the center of ICIO (a region in 
common with lPemers et al.l l l2004h ) ar e broadly cons i stent w ith the 
power law distributions discussed in IPemers et al.l l l2004ll . How¬ 
ever, we also trace these populations to even larger radius due to the 
larger field of view of MegaCam compared to CFHT12K, implying 
an even larger/more extended system as discussed above. The larger 
radial baseline allows us to show that the AGB population is more 
co ncentrated overall th an the RGB, in contrast to the conclusion 
of IPemers et al.l (l2004h . These authors also suggest that the inner 
distribution of red giants is very irregular, with a “horseshoe-like” 
shape. We are unable to comment on the very central distribution 
of red giants in our data, since the natural seeing of the site and 
the very bright starburst population in this region means we have 
severe crowding/incompleteness that affects any analysis we would 
con duct in this region. 

iNidever et al.l ( 1201 3h present strong evidence for an extended 
HI feature stretching 1.3 degrees to the northwest of ICIO. This fea¬ 
ture displays a larg e velocity gradient th at ends in a high column- 
density clump, and INidever et alJ l l2013l) argue that it is plausibly 
the result of tidal interactions in a dwarf-dwarf merger. If this pro¬ 
posed merger is also the cause of the strong star formation in ICIO, 
then the merger is relatively recent, and it is reasonable to sup¬ 
pose that the merger could leave an imprint on the stellar distribu¬ 
tion. However, no such features are observed in the optical and NIR 
datasets analysed here, and there are no obvious correlations with 
the morphology of the gas. For example, the stellar distribution is 
not obviously aligned with the direction of the extended HI, nor are 
there obvious extensions in the stellar contours that could be a sign 
of a recent interaction (either through distortion of the main body 
of IC10 or through stars left by a tidal stream from the companion). 
However, we also stress that the current optical and NIR data do not 
cover a large enough spatial extent to match with the position of the 
proposed progenitor (the optical data extends ~ 0.5 degrees away 
from the center of ICIO, whereas the HI extends for 1.3 degrees). It 
is clearly worthwhile to obtain both deeper and wider data to fully 
explore the interaction scenario and even potentially uncover a new 
companion to ICIO. Our current dat a covers a similar radial e xtent 
of earlier HI studies, for example bv lWilcots & Milled ( Il998h . and 
we refer the reader to their Figure 2 for a column density map of 
the HI in ICIO. Again, no clear correlations are present between 
the structure in the HI and the (lack of) structure in the older stellar 
distributions. 


6 SUMMARY 

In this paper, we have used optical and NIR wide field imaging 
data of the nearby. Local Group star-bursting, blue compact dwarf 
galaxy ICIO to investigate the outer stellar structure and older stel¬ 
lar populations of this unusual and highly obscured galaxy. Us- 
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ing CFHT/MegaCam and UKIRT/WFCAM we have examined the 
spatial distribution of the various stellar populations present, trac¬ 
ing them to large radius and low surface brightness. A large num¬ 
ber of bright AGB stars prevent us from obtaining a precise dis¬ 
tance to this galaxy, but we constrain it to lie in the range of 700 
- SOOkpc, consistent with earlier studies. Stellar populations of 
decreasing average age are increasingly spatially concentrated in 
ICIO, in line with results for many other Local Group galaxies and 
consistent with expectations for blue compact galaxies that have 
the star-forming region embedded within a much larger, lower sur¬ 
face brightness component of older stars. The youngest populations 
(those of the star burst and the AGB stars) are significantly offset 
from the geometric center of the older RGB population, implying 
that recent star formation occurred significantly away from the cen¬ 
ter of the galaxy. The scale of this offset is significant such that 
dynamical models of the evolution of ICIO should distinguish be¬ 
tween the center of light of this galaxy and its center of mass. The 
RGB population falls off slower than an exponential and has an 
effective radius of approximately 5.75 arcmins (~ 1.2 kpc), mak¬ 
ing this blue compact galaxy one of the most extended of the 
known dwarf galaxies in the Local Group. Importantly, the out¬ 
ermost isophotes of this galaxy are very regular in shape and es¬ 
sentially circular in morphology, and we find no evidence in the 
extended stellar component to su ggest that ICIO has u ndergone a 
recent interaction as proposed by iNidever et 'Zl j2013h . However, 
the footprints of the optical and near-IR data do not cover as large 
an area as the HI data, and it is possible that stellar signatures of 
this proposed encounter exist outside of the survey region, demon¬ 
strating the need for deeper and wider surveys of ICIO. 
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